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We use the BLM procedure to eliminate the renormalization scale am- 

O " 

biguity in the evolution equation for the non-singlet deep-inelastic structure 
function i ? 2 ,s (x, Q). The scale of the QCD coupling in the MS scheme has the 
Q^' form Q*(x) = Q(l — x) 1 / 2 jx f(x), where x is the Bjorken variable and f{x) 

cu : 

is a smoothly varying function bounded between 0.30 to 0.45. Equivalently, 
' the evolution of the n\\\ niomenl of the structure function should contain an 

effective Aqcd pattern, with A n ~ n 1 / 2 . This variation of A n agrees with 
experimental data. 
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I. INTRODUCTION 



One of the most serious problems preventing precise empirical tests of QCD is the am- 
biguity of the renormalization scale /i of perturbative predictions. Formally, any physical 
quantity should be fi independent. However, in practice, spurious /i dependence appears 
as one can only have finite-order predictions in a perturbative theory. Unless one specifies 
the argument of the coupling a s (/x) in the truncated predictions, the range of theoretical 
uncertainty can be much larger than the experimental error. Although the uncertainty of 
the prediction contributed from the renormalization scale ambiguity is expected to reduce 
as the order of the prediction gets higher, we still need a rational way of choosing a value 
for fi such that the truncated prediction will approximate the true prediction as close as 
possible. 

A conventional way of setting the renormalization scale in perturbation theory is setting 
/i to be the momentum transfer or energy scale Q of the system (provided there is only one 
energy scale); this eliminates large logarithmic terms ln(Q///) and so gives a more convergent 
series. However, often an order-1 variation of fi/Q leads to a significant uncertainty in the 
prediction of perturbative QCD. 

In multi-scale processes, the problem of scale setting becomes compounded, since the 
choice for /i can depend on any combination of the available physical scales. An example 
of this is the proper scale for the running coupling constant that appears in the DGLAP 
evolution equation for the deep-inelastic structure function. In addition to the momentum 
transfer Q of the lepton, the physical scale can also depend on the Bjorken ratio x = 
Q 2 /(2p ■ q). Equivalently, the scale controlling the evolution of each moment M n (Q) of the 
structure function can depend on both Q and n. 

Collins I], Neubert 0, and Lepage and Mackenzie have emphasized that the renor- 
malization scale should not be fixed by an ad hoc procedure but rather, should be determined 
systematically as the mean virtuality of the underlying physical subprocess. From this point 
of view, the choice of the renormalization scale \i for a particular prediction to a certain 
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order thus depends on the specific experimental measure and the truncation order. 

In this paper, we discuss how to obtain the optimal scale for the structure function 
evolution using the BLM (Brodsky-Lepage-Mackenzie) 0] scale-fixing procedure. In this 
procedure, the vacuum polarization diagrams that contribute to the non-zero QCD beta 
function are resummed into the running coupling. More technically, we absorb into the 
running coupling all factors of the number of flavors n / that appear in the coefficients of the 
perturbative expansion. This criterion automatically sets the scale to a value that reflects 
the average gluon virtuality of the subprocess. 

On the one hand, we can calculate the evolution of the moment by fixing the scale as a 
function of Q and x before integrating over x to obtain the moments; on the other hand, we 
can set the scale from the moment-evolution equation for each n. We give the appropriate 
scale for the moment-evolution equation for each n and show that the two procedures give 
the same result to the order that we consider. 

II. THE BLM SCALE-FIXING PROCEDURE 

The crucial idea of BLM procedure {||| is that, in the ay scheme [ay is defined through 
the heavy quark potential V{Q 2 ) by V{Q 2 ) = —4:7iCf(Xv(Q)/Q 2 ], the scale must be chosen 
to absorb all the vacuum polarization (non-zero beta function) contributions into the run- 
ning coupling constant. Thus one has to choose the argument of the coupling constant in 
each order of perturbation theory so that there is no n/ dependence in the coefficients of the 
coupling constants (those light-by-light contributions that are not associated with renormal- 
ization are not resummed). This is a good choice of scale since all the vacuum-polarization 
contributions are then automatically summed to all orders for any finite-order prediction. 
For a next-to-leading order (NLO) perturbative prediction in any renormalization scheme, 
the procedure simply translates to choosing the scale such that there is no nj dependence 
in the NLO coefficient of the coupling constant, by the transitivity of the procedure to this 
order [[5]. 
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To show explicitly how to set the scale for an NLO perturbative QCD prediction, consider 
a prediction in any scheme of the following form, 



P = Po<x<Ji) 1 + [A(fi)n f + B(ji)\ 



a(p) 



+ 0{or 



71 



(1) 



where p is the renormalization scale. The parameters p, po, A and B are all in general 
dependent on one or more physical scales such as center-of-mass energy, momentum transfer 
and Bjorken parameter x. The p dependence of the coefficients A(p) and B(p), fixed by the 
renormalization group 

dp 



dp 



are of the form 



+ O(a d 



1 



(2) 



A(p) = A'-- ln(p) 

B(p) = B' + ^C A Hp) , 
o 



(3) 



where A' and B' are p independent. 

We apply the BLM procedure to eliminate the rif and /3 dependence in Eq. ([I]). We 
choose the renormalization scale p to be 



Q* = pexp(3A(p)), 



(4) 



so that 



A{Q*) = 

[note that the right-hand side of Eq. (Q) is p independent by Eq. (|3])]. 
Using the BLM scale, the perturbative QCD prediction becomes 



(5) 



P = Poa(Q* 



l + B(Q* 



7T 



+ 0(a d 



(6) 



Rewrite it in terms of the original coefficients A(p) and B(p) using Eq. 



P = Poa(Q*) { 1 + 



^C A A(p) + B(p) 



a(Q* 
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+ 0(a 3 ). 



(7) 



This is the scale- fixed perturbative QCD prediction. With the chosen scale Q*, all 
vacuum polarization is resummed into the running coupling constant, with the coefficient of 
a 2 independent of nj. 

In short, by comparing Eq. (|IJ) and Eq. (|7|), setting the scale for an observable is simply 
equivalent to replacing rif by yC^, and using Q* as the argument of the coupling constant 
(/i dependences in A(fi) and B(/i) cancel). 

For a prediction with \i preset to any energy scale Q, as is usually given in the literature, 
the procedure for resetting the scale appropriately is identical to the one presented above 
except with Q replacing \i. 

As we have seen, to obtain the renormalization scale only the rif term of the NLO result 
is needed. Therefore, one can improve a leading-order result by setting an appropriate scale 
without calculating the full NLO prediction. Extensions of the method to higher orders are 
given in [|, @, and 0. 

III. SCALE SETTING OF THE STRUCTURE-FUNCTION EVOLUTION 

EQUATION 



The non-singlet, structure- function evolution equation is 



8T m 



x y 



x 



x 



P qq [-,a(Q)\±Pg q [-,a(Q) 



(8) 



where, in terms of experimental quantities, 
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(Ff - Ff)/x for+ , 
(Ff - Ff)/x for- , 

and P qq ± Pqq are the non-singlet evolution kernels, for which the NLO perturbative predic- 
tions in MS scheme with the renormalization scale preset to momentum transfer Q are |§, 
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+6(1 -x) ( dxPg q (x, a m (Q)) + 0(a 3 s ). 
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P, q (x,a m (Q))= (^f^) (Cf -~C a )C f [2(1 + x)lnx + 4(1 - x) 



with 
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(10a) 



(10b) 



T = n f /2 

C A = 3 

CV = 4/3 

Li 2 (x) = - r^ln(l-y). 
Jo y 

The + distribution notation serves as a regulator defined as 

#(+ 



dxf(x) 



1 — x 



JO 



(11) 



(12) 



1 — x Jo 1 — X 

Since all the coupling constants reside in the evolution kernels, we set the scale for the 
structure-function evolution equation by setting the scale for the kernels. When the kernels 
P qq ± Pq q have been recast in the form of Eq. ([J), the coefficient of rif can be obtained as 

l-x\ 29 (l-x 



A(x) 



In 



x- 
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1 13 

1 + x 2 J V4 + ~4 X 



(13) 



Using Eq. (f4|), BLM scale- fixing procedure gives for the BLM scale 

Q*(x) = Qexp(3A(x)) = Q (1 - ^ f(x), 

x 



(14) 



where 



f(x) = exp 



l-x\fl 13 \ 29 



(15) 



is a smoothly varying function bounded between 0.30 and 0.45 as shown in Fig. |I[ 

Note that the above manipulation is done inside the + distribution notation. It is because 
after setting the scale to be dependent on x, the coupling constants must also be included 
in the + distribution notation in order to preserve the Adler sum rule. 

Finally, the scale-fixed evolution kernels can be obtained by rewriting Eq. ([H]) in the 



form of Eq. (0), using the BLM scale given by Eq. ([[J) 



P„(x,a m (Q'( X ))) = < "MsW-M^ fl+f^ + MW 
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1-x J 2 V ; 

i±f!) L* x+ m_«l\ 
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+ 2(l + x)lnx-6x + 6]}}_ 



+ 5(1 - x) / dxP qq (x, a m (Q*(x))) + 0(ai), 

J 



(16a) 



and 



P qq (x,a m (Q*(x)))= ^ ^MsgM ^j (C F -ic A )C F [2(l+x)lnx + 4(l-x) 

1 +X 2 / 7T 2 ~ 

H In 2 x - 41nxln(l + x) - 4Li 2 (-x) 

1 + x V 3 



(16b) 



Q 



The scale fixed-structure function evolution equation is then given by 
dF? s , _ f 1 dy 



2 W ' / 2 



(x,Q) 



x y 



^ s (y,Q). (17) 



This equation can also be written in terms of ay using the relation 



a Ms(^) = ov(/-*exp(5/6))(l + 2a v /ir + ...). 



(18) 
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IV. VERIFICATION OF THE SCALE 



A good way to see the validity of the scale set for the structure-function evolution equa- 
tion is to apply it to moment evolution and compare the results with experimental data. 
The moment analysis of the Fermilab muon deep-inelastic scattering data |TT| showed an 
interesting variation of the effective A QCD with n, both at leading and next-to-leading order. 
Here we would like to demonstrate that this variation can be readily explained by the choice 
of the renormalization scale even at leading order. 

The nth ordinary moment of the non-singlet structure function is defined as: 

M„(Q 2 ) = f 1 dxx n ~ l ^ s (x, Q 2 ). (19) 

The Nachtman moments will give the same result up to higher twist terms. 

Taking the nth moment of the structure function evolution equation (Eq. (|17D ), inter- 
changing the order of integration on the right-hand side and rearranging, one can obtain the 
equation of evolution for the moment: 

d \*t<§ 2) = L dxXn ~ l { P m{^oc{Q*{x)))±P m {x,a{Q\x))) , (20) 

which is just the nth moment of the evolution kernels. 

Let us consider only the leading order prediction of the kernel function in Eq. (|16a| ), 
then Eq. ( |20"D becomes 

here Q*{x) is given by Eq. ([TJ]), and 

P» = Cf 1 -^ (22) 

the leading order coefficient. 

The crucial difference between the conventional calculation, that is, using the momentum 
transfer Q as the renormalization scale, and the calculation here, namely using the BLM 
scale, is that the argument of the coupling constant in Eq. (|2lD is now x dependent. 
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To see the effect of the scale setting on the leading-order analysis, expand a(Q*) in 
powers of (3 a(Q) (similar results can be obtained without the expansion), that is, 



a(Q*(x)) = a(Q) {l + fa In ( J^) + . . . 



then Eq. Q2"I|) becomes 

d In M n (Q 2 ) a(Q) 



<91nQ 2 



2vr 



-4.„ 



Q*(x) J 2tt 
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(23) 



(24) 



where 



.4 



„- , dxx 11 1 P g ° g (x)+,and 



B n = ^- f 1 fa"" 1 
A n Jo 



(25) 
(26) 



Now Eq. ( |2"4D can be integrated with respect to InQ 2 by using, to LO (to compare with 
LO experimental analysis), 

47T 



a(Q) 



/3 ln(g 2 /A 2 ) 



(27) 



and the solution for M n (Q) is then 

M n (Q)~ 1/dn = C n ln((5 2 /A 2 )exp 



2B r . 



ln(g 2 /A 2 ) 



(28) 



where d n = —2A n / (3 Q and C n are constants. 

This is our prediction for the evolution of the moments. 

For experimental data taken at large Q, the exponential term can be expanded and we 
have an equation linear in ln(Q 2 /A 2 ): 



M n {QY 1 ' d - = C ri [ln(Q 2 /A 2 ) - 2B 



Had we set the scale to Q, we would have had 



(29) 



M^Q)- 1 '^ = C n [ln(Q 2 /A 



(30) 



where A n is an "effective" A at LO and will depend on n as predicted by Eq. (p9|) as 



A n = A/exp(-S n ). (31) 

Values of exp(— B n ) for n from 2 to 10 are plotted in Fig. |2|; these are well-described by 
the form 

exp(-5 n ) = 0.656/n 1/2 . (32) 
Therefore, we expect the dependence of A n on n to be 

n l/2 

A » = a«6 A - < 33 > 

In fact, the data of Gordon et al does show the predicted variation. In Fig. ^ we show 
the results extracted for A n against n. A fit to Eq. ( p3|) yields A^jg = 170 ± 120 MeV, where 
the error is estimated from higher-order terms in Eq. 



V. SCALE SETTING FOR THE MOMENT-EVOLUTION EQUATION 

The BLM scale fixing procedure can also be applied to the moment-evolution equation 
after the moment integration has been carried out, with evolution kernels given by Eq. (P~D|) . 
When the resulting equation is written in the form of Eq. (|l|), the coefficient of nj is given 
by 

A = -p 1±, (34) 



where A(x) is given by Eq. (13). This is related to our previously defined B n as 



A = (35) 

Again, using Eq. (f4|), the BLM scale for the moment-evolution equation is 

Q* n = Qexp(-B n ). (36) 

Now the scale is x independent but n dependent. 

The scale for the moment-evolution equation can also be read from Fig. ^| and has a fit 
of Qn/Q — 0.656/n 1 / 2 . The n 1//2 behavior for large n was also predicted in [|J. 
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Keeping only the leading-order term of P qq , the moment-evolution equation can be writ 
ten analogous to Eq. (^|) as 

d\nM n (Q 2 ) a(Q* n 



<91nQ 2 2tt ~ An ' 

Again, it can be solved by rewriting the coupling constant using Eq. (^7|) to be 

M n {Q)~ l ' d - = C n \n(Q* n 2 /A 2 ) 



(37) 



\n(Q 2 /A 2 ) - 2B n 



(38) 



which is simply Eq. (^). Actually, the commutativity of scale setting and moment integra- 
tion works at any higher order, because a choice of the renormalization scale can only affect 
the result by an order higher than that of the calculation. 



VI. CONCLUSION 

We have applied the BLM procedure to set the scale in the evolution equations for the 
non-singlet structure functions and their moments in deep-inelastic scattering. The variation 
of the effective A obtained from the nth-order moment-evolution data is explained by this 
method, and an unambiguous value of Aj^g is obtained. 
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Fig. 2 

FIG. 2. Plot of exp(— B n ) (or Q^/Q) against n. The fitted curve has the form 
exp(-B n ) = 0.656/n 1 / 2 . 
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Fig. 3 

FIG. 3. Plot of A n against n, after B. A. Gordon et al. 
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